Phonons in Superconducting CaCe Studied via Inelastic X-ray Scattering 
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We investigate the dispersion and temperature dependence of a number of phonons in the recently 
discovered superconductor CaCg utihzing inelastic x-ray scattering. Four [OOL] and two a6-plane 
phonon modes are observed, and measured at temperatures both above and below Tc. In general, 
our measurements of phonon dispersions are in good agreement with existing theoretical calculations 
of the phonon dispersion. This is significant in light of several discrepancies between experimental 
measurements of phonon- derived quantities and theoretical calculations. The present work suggests 
that the origin of these discrepancies lies in the understanding of the electron-phonon coupling in 
this material, rather than in the phonons themselves. 



PACS numbers: 74.25.Kc,74.70.-b,78.70.Ck 

CaCe has attracted a great deal of interest since the 
recent discovery of its superconductivity below 11.4 K 
P, The transition temperature is unprecedented in 
graphite intercalation compounds, where transition tem- 
peratures are, with one exception, an order of magnitude 
lower. (The second highest known transition tempera- 
ture, of YbCe, is still nearly 50% lower, 6.5 K.) Despite 
an intriguing initial suggestion of a novel pairing mech- 
anism [3| much of the early experimental experimental 
evidence indicates that CaCe is a conventional weakly- 
coupled superconductor [4, 5]. Still, the unusually high 
transition temperature has prompted a great deal of spec- 
ulation as to which phonons play a part inpairing and 

why the transition temperature is so high iiiaii. 

Calandra and Mauri Q and Kim et al. [4] have cal- 
culated the phonon spectra and the electron-phonon 
coupling, using density functional theory (DFT). Both 
predict electron-phonon coupling with the Cac (out-of- 
plane) and Cab (in-plane) modes. While these two calcu- 
lations of the phonon dispersions are similar, their pre- 
dictions of the coupling differ, with Calandra and Mauri 
assigning additional coupling to the Cab modes. 

There are a number of discrepancies between the theo- 
retically calculated phonon spectra and electron-phonon 
coupling and a variety of phonon-derived experimental 
quantities for this system - as outlined in a recent paper 
by Mazin et al. [TO]- For example, recent measurements 
of the isotope effect, a{X) = —{d\ogTc/dMx) (where Tc 
is the superconducting transition temperature and Mx 
is the mass of atom X) have found a(Ca) = 0.53 [ll|, in 
disagreement with the calculated value of a(Ca) = 0.24 
jsj. This suggests that the coupling to soft Ca modes is 
greater than predicted by calculation. Furthermore, the 
measured temperature dependence of the specific heat 
does not show the large coupling effects derived from 
the calculated electron-phonon spectral function, a'^F{uj) 



0, [T^. (Also called the Eliashberg function, a'^F{u;) 
is a measure of the electron- pho non coupling. F{u;) is 
the phonon density of states [i2[-) Thus, in contrast to 
the isotope effect measurements discussed above, specific 
heat experiments suggest that the coupling to acoustic 
Ca modes is significantly weaker than is predicted by 
theory. 

Additionally, the upper critical field of CaCe has been 
shown to have a linear temperature dependence down to 
1 K while calculations of the upper critical field, based on 
LDA calculations of have predicted non-linear 

behavior below 4 K BElil- Finally, DFT calculations 
of the phonon spectra understate the increase of Tc under 
pressure by an order of magnitude [4, 14, 15]. 

The measurements of the isotope effect and the specific 
heat appear to lead to contradictory conclusions, since 
the isotope measurement implies larger-than-calculated 
coupling to soft Ca modes while the specific heat mea- 
surements suggests less-than-calculated coupling to soft 
Ca modes. This apparent conflict rests on a comparison 
of experimental results to theoretical predictions which 
yield A - 0.85 0, 0. However, a recent measurement 
suggests that the superconductin g ga p is 2.3 meV, 40% 
larger than previously measured [16]. Such a large gap 
suggests that CaCe is a strongly coupled superconductor. 

It is clear that the theoretical understanding of CaCe 
is by no means complete, and in particular a number 
of experimental observations of quantities that are ei- 
ther directly or indirectly related to the phonon spectra 
are inconsistent with the present theoretical description. 
Given this situation it is clearly of great interest to mea- 
sure the phonon spectra and compare it to the existing 
theoretical predictions [4, 6] . 

In this paper we report measurements of the phonon 
dispersions in CaCe, utilizing inelastic x-ray scattering 
(IXS) to do so. In particular, we have measured the dis- 



persion, and temperature dependence of the dispersion, 
for a number of the low energy modes, including soft Ca^b 
modes and part of the Cab mode, that are believed to be 
important in the superconductivity of the CaCe [6, 7]. In 
general, we find good agreement with the extant predic- 
tions for the phonon dispersion curves. These results ap- 
parently confirm the validity of the DFT approach used 
in such calculations and suggest that the origin of the var- 
ious discrepancies with experiment discussed above lie in 
the understanding of the electron-phonon coupling not of 
the phonon modes themselves. 

The samples were prepared via immersion of a Highly 
Orientated Pyrolytic Graphite (HOPG) platelet in a 
lithium/calcium alloy for 10 days, as described in detail 
by Pruvost et al. (17]. The HOPG was ZYA grade pur- 
chased from GE Advanced Ceramics with an initial (non- 
intercalated) highly- aligned c-axis mosaic of ~ 0.4°. It 
is powder-like in the a6-plane. The lithium and calcium 
were purchased from Sigma- Aldrich at 99.99% purity. X- 
ray diffraction of the resulting shiny silver CaCe platelet 
is shown in Fig. [it. Diffraction showed very high sam- 
ple purity with the intensity of the (002) Bragg peaks of 
the LiCe phase and the graphite phase having less than 
0.25% the intensity of the CaCe (006) Bragg peak. 

The sample dimensions were 4 x 4 x 0.7mm with a 
post-intercalation c-axis mosaic of ~ 3.5°. The crystal 
structure and Brillouin zone of CaCe (space group i?3m, 
a=5.17 A, a = 49.55) are shown in Figs.[Tt and b [2]. 

After synthesis, the samples were mounted in a beryl- 
lium dome in an argon or helium atmosphere. The region 
around the beryllium dome was subsequently pumped 
out to rough vacuum. The diffraction pattern of the sam- 
ple was checked before, after and several times during the 
experiment. Both the out-of-plane diffraction patterns 
and the rocking curves of the (003) and (006) refiections 
were monitored during the experiments and did not show 
any signs of degradation [l^] • A post-experiment diffrac- 
tion pattern is shown in Fig. [it. 

The IXS experiments were performed at sector 3 at 
the Advanced Photon Source of Argonne National Lab- 
oratory [2O, [lij. Data were collected using Si(18 6 3) 
analyzer reflections. The instrument had an overall reso- 
lution of 2.3 to 2.5 meV. Photon flux was 1.8 - 2.4 x 10^ 
photons/second. Four analyzer crystals and four inde- 
pendent detectors allowed data collection at four mo- 
mentum transfers simultaneously. The momentum res- 
olution was 0.072 A in the scattering plane and 0.18 
A perpendicular to it. The beam size at the sample was 
250 X 300/im. The spectra were normalized by a beam in- 
tensity monitor immediately before the sample. Spectra 
with Q in the [OOL] direction were measured in refiection 
geometry while spectra with Q in the a6-plane were mea- 
sured in transmission geometry. Finally, we note that 
identical results were obtained from two samples, pre- 
pared at different times and measured in separate runs. 

A series of spectra taken at 5K with Q in the [OOL] 
direction and in the fourth Brillouin zone, with Q be- 



(a) 





(c) 



CaCfi 



^^^^^ LiCg Graphite 
i (002) (002) 



CaCg 
(006) 



LiCe 
(004) 



1.5 



2.0 2.5 
Q [OOL] {A-1) 



3.0 



FIG. 1: (color online) (a) The CaCe crystal structure 0- 
Green circles represent Ca atoms and the grey ring is graphene 
jJi]. (b) The CaCe Brillouin zone, (c) An [OOL] diffraction 
pattern taken after the inelastic experiment. The pattern 
exhibits sharp CaCe (003) and (006) peaks. The locations 
of possible impurity peaks from LiCe and Graphite are indi- 
cated. Almost no scattering is observed at these locations. 



tween |G and 2G where G is one reciprocal lattice vector, 
are shown in Fig. [2l Typical phonon intensities are < 1 
count /second. Example fits are shown in the figure. The 
elastic line is fit with a psuedo-voigt while the phonon 
peaks are fit with Lorentzians. Near the L point zone 
boundary, the feature at 10 meV is seen to broaden and 
a better fit is obtained if the data are fit to two modes. 

A summary plot showing the resulting peak positions 
is shown in Fig. [3l Note that not all of the phonon ob- 
servations reported in Fig. [3] are visible in Fig. [2] because 
the cross-section of different phonons varies with momen- 
tum transfer. Specifically, the phonon intensity is propor- 
tional to {Q ' e)'^ / hco ^ where e*is the atomic displacement, 
Q is the momentum transfer and hoo is the energy of the 
phonon mode. Some phonons shown in Fig. [3] are visi- 
ble only in higher Brillouin zones or at longer counting 
times. The data in Fig. [3^ are shown in a reduced zone 
scheme with spectra from different Brillouin zones folded 
back into the first Brillouin zone. 

We discuss first the F — L dispersion (Fig. [3^). Four 
phonon bands are observed in this direction. The low- 
est energy band, labelled A, is almost dispersionless over 
the measured region and has an energy near 10 meV. 
The second band, labelled B, disperses from meV to 
25 meV. A small portion of a third band, labelled C, can 
be seen near L near 10 meV. The final observed band, 
labelled D, disperses between 25 meV and 40 meV. 

It is difficult to assign an exact cross-section to a par- 
ticular mode because the data were collected in differ- 
ent geometries. The elastic intensity is a reasonable, 
but imperfect, measure of the amount of material in 
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FIG. 2: (color online) Inelastic spectra along F — L from the 
4tli Brillouin zone. As discussed in the text, the phonon cross- 
section means that some of the phonon modes reported in 
Fig. [3] do not appear in these data. Data are labelled by 
momentum transfer in the reduced zone scheme. The spectra 
are offset for clarity and are normalized to account for the 
different efficiencies of the different detectors. The y-axis is 
a log^Q scale. The blue (dark gray) long dashed lines are 
the total fit. Each red (light gray) short-long dashed line is 
a portion of the total fit, including the elastic line and one 
phonon mode peak. 



the scattering volume, and additional factors complicate 
the comparison of data taken in the reflection ([OOL]) 
and transmission (a6-plane) geometries. However, we 
can report the relative cross-sections of various modes 
at a particular momentum transfer in the [OOL] direc- 
tion. At Q =2.08 A , which corresponds to a reduced 
Q =0.69 A , modes C and D have roughly equal intensi- 
ties of 0.5 counts/second while mode A has an intensity 
of roughly 0.1 counts/second at 5K. These count rates 
represent the raw data and are not scaled by the IXS 
polarization factor, (Q • ef' jTiuj. 

Also shown in Fig. [3^ are calculations of the phonon 
bands [6]. The two calculations agree to within a few 
meV in the measured energy range with the exception of 
the dispersion of the lowest energy optical band in the 
a&-plane 0, @]- The calculation of Calandra and Mauri 
^ agrees with the data somewhat better than the cal- 
culation of Kim et al. [4], so it is shown here. Moreover, 
Calandra and Mauri identifj^the phonon mode associated 
with each calculated band [6]. By comparison of our ex- 
perimental data with this calculation, we identify mode 
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FIG. 3: (color online) Phonon dispersion from spectra taken 
at three different temperatures, 300 K in black circles, 05 50 
K in red triangles. A, and 5 K shown in blue squares, □. (a) 
shows the [OOL] spectra. Calculated dispersions are shown as 
dashed lines [y]. Different calculated modes are labeled by 
A-E. (b) shows the a6-plane spectra. The distance to x is 
indicated by a vertical black long dashed line. The F-X cal- 
culated dispersions are shown as pink short-long dashed lines 
and the F-x dispersions are shown as blue dashed lines. Note 
that F-X is not a strictly a5-plane dispersion. The theoretical 
curves most closely corresponding to the data shown marked 
with thicker lines than other curves. 



A as a primarily G^ah oscillation with some Cc character, 
mode B as a mixed Cac and Cc mode, and mode D as a 
mixed Cac and Cc mode [6]. 

In general, the agreement between theory and exper- 
iment is good, though we note that the calculations of 
modes A and D understate the energy of the measured 
band throughout the dispersion. In addition, there are 
predicted modes which are largely or entirely absent in 
our data: modes C and E. This apparent absence does 
not necessarily mean that the mode does not occur. In 
particular, mode C is predicted to have primarily Ca^b 
character, which, because of the phonon polarization fac- 
tor in the IXS cross-section, means it should not be ob- 
served. This mode also has a small amount of Cc charac- 
ter, in addition to the majority Ca^b component. A small 
Cc amplitude may be why these measurements observe 
the C mode in only a small portion of reciprocal space. 
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In fact, mode A is also predicted to have primarily Caab 
character and a small amount of Cc character, consistent 
with this. Near Q =1.24 A the measured intensity of 
this mode is 45 times larger in the a6-plane than it is in 
the [OOL] direction (after scaling by the elastic intensity) . 
Thus, the polarization factor also explains the asymme- 
try in the intensity of this mode. 

Despite measuring a wide range around 50 meV we are 
not able to observe 50 meV Cc phonons modes. Spectra 
were counted for between 80 and 120 seconds per 0.25 
meV step. We believe that a 10 count peak would have 
been observable, which places an upper bound on any 
such mode of 0.1 counts/second at 300K, less than 10% of 
the dispersing phonon band at similar momentum trans- 
fers. A related Cc phonon mode was also not observed in 
an IXS measurement on graphite Q . We speculate that 
the amplitude of these modes is very small or the lifetime 
is very short (resulting in a broad peak). Either of these 
effects, when combined with the high predicted energy of 
mode E, would result in a very small IXS intensity. 

Next we turn to the a6-plane phonons (Fig. [8)3). The 
interpretation of the a6-plane phonons is somewhat com- 
plicated because the data are not taken at a point in 
reciprocal space. Rather, because the sample is not ori- 
ented in-plane, each spectra is taken at fixed |Q|, which 
represents an average over a circle in the a6-plane in re- 
ciprocal space of radius |Q|, centered at P. Therefore, 
phonons which have a significant anisotropic dispersion in 
the a6-plane are not visible in these experiments because 
this anisotropy will translate into an energy broadening 
of the measured mode and consequently a reduction in 
the peak intensity. Of course, if the a6-plane dispersion is 
isotropic (identical in all a6-plane directions) the average 
is irrelevant and the mode may be observed. Finally, all 
the data must be taken in the first Brillouin zone because 
the powder structure of the a6-plane and the unequal 
distances of different a6-plane directions in the Brillouin 
zone mean that phonon positions cannot be folded back 
to the first Brillouin zone. Despite these difficulties two 
phonon bands are observed. The flat band near 10 meV 
and a highly dispersive band which we measure between 
15 meV and 50 meV. At Q =0. 58 A" ^ and 5K, the two 
measured a6-plane phonon modes have roughly equal in- 
tensities of 0.1 counts/ second. 

Comparison with Calandra and Mauri (Fig. [3]) iden- 
tifies the low energy band as a principally Ca^b mode 
and the dispersing band as a Cab mode [6|. This latter 
mode, which is the a6-plane extension of mode C mea- 
sured in the [OOL] direction, discussed above, disperses 
to higher energy than predicted [4] . A second highly dis- 
persing Cab mode is predicted by Calandra and Mauri, 



but is not observed . This may be because the band is 
not isotropic. In fact, the principal difference (below 60 
meV) of the calculations of Calandra and Mauri and Kim 
et al. is that Calandra and Mauri predict two isotropic 
Cab bands dispersing from F at 15 meV while Kim et 
al. predicts that only one of them is isotropic. 

Calandra and Mauri [6] and Kim et al. [3] both pre- 
dict a flat, isotropic band near 38 meV. This mode is 
not observed in experiments, which cannot be due to the 
a6-plane averaging effects discussed above. Thus, its ab- 
sence in the experimental data is explained by its assign- 
ment as a mixed Cac and Cc mode: the IXS polarization 
factor would suppress any such mode in this geometry. 

Finally, the phonon dispersion is measured at three 
temperatures 5 K, which is below the 11.4 K transi- 
tion temperature; 50 K and 300 K. No differences in the 
phonon spectra at different temperatures are observed, 
within the limits of our measurements. Additionally, 
there are no significant changes in the widths. Also, the 
intensity of the peaks appears to follow the temperature 
dependence predicted by Bose statistics, as expected. 

In summary, we find good agreement between exper- 
iment and theory, although calculation underestimates 
the energy of two modes in the [OOL] direction. This 
agreement has a number of implications for the gen- 
eral understanding of the system. One immediate conse- 
quence is that the predicted F{uj), the phonon density of 
states, appears to be largely correct 0, @]- 

Nevertheless, as discussed in the introduction, serious 
disagreements exist between experiment and the theo- 
retical understanding of the system. The slight disparity 
between the calculated and measured F{u;) is probably 
not sufficient to reconcile these results. The disagree- 
ments are largely a result of predictions made using the 
calculated a'^F{uj) 0, Q. Thus, by confirming our un- 
derstanding of the phonon dispersion, and hence the cal- 
culated F{uj)^ the present work raises serious questions 
about the calculation of the electron-phonon coupling. 
We hope that this work inspires theoretical studies to rec- 
oncile the existing discrepancies with experimental data. 
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